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^ : Abstract 



We have performed a critical analysis of different approximation schemes 
for the calculation of two-neutrino double beta decay (2/32z/-decay) matrix 
elements. For that purpose the time integral representation of the 2(32u- 

X' 

' decay matrix element has been used. We have shown that within the single- 

particle approximation of the nuclear Hamiltonian the 2/32z/-decay matrix 
element is equal to zero because of the mutual cancelation of the direct and 
cross terms. The quasiboson approximation (QBA) and renormalized QBA 
(RQBA) schemes imply for the 2/32i/-decay transition operator to be a con- 
stant, if one requires the equivalence of initial and final Quasiparticle Random 
Phase Approximation (QRPA) and renormalized QRPA (RQRPA) Hamilto- 
nians. It means that 2/52i^-decay is a higher order process in the boson expan- 
sion of the nuclear Hamiltonian and its higher order boson approximations 
are important. The equivalence of the initial and final QRPA and RQRPA 



Hamiltonians has been discussed within the QBA and RQBA, respectively. 
We have found that the mismatching of both Hamiltonians is getting worse 
with increasing strength of particle-particle interaction especially in the case 
of QRPA Hamiltonians. It is supposed to be one of the reasons of the ex- 
treme sensitivity of the 2/?2z/-decay matrix element to the residual interaction 
appearing in explicit calculations involving the intermediate nucleus. Fur- 
ther, the Operator Expansion Method (OEM) has been reconsidered and new 
2/32z/-decay transition operators have been rederived in a consistent way. The 
validity of the OEM approximation has been discussed in respect to the other 
approximation schemes. The OEM combined with QRPA or RQRPA ground 
state wave functions reflects sensitively the instabilities incorporated in the 
considered ground states. Therefore, the predicting power of the OEM should 
be studied with help of the other ground state wave functions, e.g. shell model 

ones. 
23.40.Hc 
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Typeset using REVT^ 



I. INTRODUCTION 



The two-neutrino double beta decay (2/32z/-decay) is a second order process of the weak 
interaction, which is allowed by the Standard model [0]. In 2/32z/-decay the nucleus (A,Z) 
undergoes the transition to nucleus (A,Z+2) with emission of two electrons and two antineu- 
trinos. This rare process is already well established experimentally for a couple of isotopes. 
The inverse half-live of 2/32i/-decay is free of unknown parameters on the particle physics side 
and is expressed as a product of a phase-space factor and the relevant 2/52z/-decay nuclear 
matrix element. Since the phase-space factor can be calculated with the desired accuracy, 
the experimental studies of 2/52i/-decay give us directly the value of the 2/32z/-decay nuclear 
matrix element. In this way 2/?2z/-decay offers a sensitive test of nuclear structure calcula- 
tions. The calculation of the 2i//5/3-decay nuclear transition continues to be challenging in 
view of the smallness of the predicted nuclear matrix elements and the fact that the mech- 
anism which is leading to the suppression of these matrix elements is still not completely 
understood. 

The proton-neutron Quasiparticle Random Phase Approximation (pn- QRPA) has been 
the most frequently used nuclear structure method for evaluating 2/?2z/-decay rates because 
of the remarkable success achieved in revealing the suppression mechanism of 2/?2i/-decay 
matrix elements However, the extreme sensitivity of the 2/32z/-decay matrix elements 

on the pn 1+ particle-particle matrix element as well as the collapse of the QRPA solution in 
the physically acceptable region of the particle-particle strength of the nuclear Hamiltonian 
renters it difficult to make definite rate predictions. 

Some attempts have been done to overcome the above drawbacks, e.g. by including 
higher order RPA corrections [^, particle number projection ||^J^ and proton-neutron pair- 
ing 1^ in the theory. However, none of these modifications of the QRPA prevents the 
collapse and inhibit the nuclear matrix element to cross zero close to the physical value of 
the particle-particle force. Recently, Toivanen and Suhonen have proposed a proton-neutron 
renormalized QRPA (pn-RQRPA) p, which goes beyond the QRPA and takes into account 



3 



the Pauli exclusion principle in an approximate way. It has been shown that the above 
phenomena could be connected with the limitation of the QRPA approach - the quasiboson 
approximation, which violates the Pauli exclusion principle. The renormalized quasibo- 
son approximation (QBA), on which the pn-RQRPA is based, inhibits the collapse of the 
pn-RQRPA solution for a physical value of the particle-particle interaction strength. In ad- 
dition, the calculated 2/32i/-decay nuclear matrix elements via pn-RQRPA have been found 
significantly less sensitive to the particle-particle force within its physical values in respect 
to those obtained by the pn-QRPA 0~ |jl2|. This behavior has been confirmed also within 



the renormalized QRPA with proton-neutron pairing (full-RQRPA) ||10|| . 

In spite of the advantages, of the renormalized QRPA over the QRPA, the RQRPA can 
not be considered as the ultimate solution for the calculation of the 2/32i/-decay process. 
Several shortcomings still plague the RQRPA: i) The Ikeda sum rule is violated ii) 
There are two sets of the intermediate nuclear states in the calculation generated respectively 
from initial and final nuclei, which do not coincide with each other, iii) The Pauli exclusion 
principle is taken into account only in an approximate way. All these leaps of faith of the 
RQRPA approach have common origin. It is the particle number non-conservation. The 
effect of these shortcomings on the 2/52i/-decay amplitude is well understood. 

A longstanding problem of large discrepancies between the values of the predicted and 
calculated 2/52z/-decay matrix elements has led to a development of alternative methods, e.g. 
one of them is the Operator Expansion Method (OEM) p3[. The OEM tries to avoid 



the necessity of evaluating the sum over the intermediate nuclear states. The price paid for it 
is that one has now to deal with more than two-body operators and commutators involving 
kinetic energy terms in the commutator expansion of the 2/?2z/-decay transition operator. 
Recently, this method has been reconsidered ||23[| and it has been shown that the previ- 
ous derivation I^Tf of the OEM-potential was not consistent. The OEM-calculation 



with a consistent OEM-potential combined with the pn-RQRPA ground state way functions 
(OEM-I-RQRPA) has exhibited a large sensitivity on the strength of the particle-particle 
force within its physical values |2^. There is a speculation that the approximations of the 



pn-RQRPA method are responsible for this behavior. Therefore, detail OEM+RQRPA cal- 
culations are expected to be helpful for solving the problem of stability of the 2/32i/-decay 
matrix element in respect to model parameters. 

Till now no consistent many body approach is available for the calculation of the many- 
body Green function governing the 2/52i/-decay process, because of the computational com- 
plexity of the problem. Therefore we can not avoid the introduction of different approxima- 
tion schemes in the evaluation of the nuclear matrix elements. Nevertheless, we can try to 
understand the limitations of the different approximation schemes. The aim of this work is 
to perform a critical analysis of the QBA and renormalized QBA schemes by using the time 
integral representation of the 2/52z/-decay nuclear matrix element and to discuss the validity 
of the QRPA, RQRPA and OEM+RQRPA calculations. 



II. 2/?2zy-DECAY NUCLEAR MATRIX ELEMENT. 

If the two-nucleon mechanism for the 2/?2z/-decay process is considered, then for the 
matrix element of this process we have 



<f\S^'^\^> 
xm(pi)7m(1 + l5)ui-ki)u{p2)'J^{l + 75)m(-^2) 

-{Pi ^ P2) - {h ^ + {pi ^ P2){ki ^ k2), (1) 

where 

J^.iPuP2,ki,k2) = J e-^(fi+fciKe-*(P2+'=2):.. 

out <Pf\T{Jf,{xi)J^{x2))\Pi >in dxidx2. (2) 

Here, Np = (l/(27r)3/2)(l/(2po)^/^), Pi and p2 {ki and ^2) are four-momenta of electrons 
(antineutrinos), pi and pf are four-momenta of the initial and final nucleus. J^{x) is the 
weak charged nuclear hadron current in the Heisenberg representation [0,^ . 



The matrix element in Eq. contains the contributions from two subsequent nuclear 
beta decay processes and 2/52i/-decay |^. They could be separated, if we write the T- 
product of the two hadron currents as follows: 

T{Jf,{Xi)Ju{x2)) = Ji,{Xi)Jy{x2) + 0(X2O - Xiq)[Ju{x2), J fji{Xl)\. (3) 

The first term in the r.h.s. of Eq. is associated with two subsequent nuclear beta decay 
processes, which are energetically forbidden for the most of 2/52z/-decay isotopes. The second 
term corresponds to 2/52i/-decay process. We see that the 2/32z/-decay nuclear matrix element 
is given by the non-equal-time commutator of the two hadron currents. It will be shown later 
that this feature is crucial for our understanding of the approximation schemes of different 
nuclear models. 

We further assume the following standard approximations: i) The non-relativistic impulse 
approximation for the hadronic current Ji/(0, y). ii) We keep only the contribution from the 
axial current, iii) Only the si/2 wave states of the emitted electrons are considered, iv) Our 
interest will be restricted only to the most favored ^tniuai ^ ^%nai nuclear transition. Then 
we have. 



vk 



x2n6{Ef- Ei+pio + ho+P2o + k2o), k = 1,2,3, (4) 



where, 



Mgt = - (e^*(Pio+'=io-A) + e't(p^o+k2o~A)^MAA{t)dt, (5) 

^ J 



with 



MAA{t) =f< 0+|^[A,(V2), Afc(-t/2)]|0+ >„ Akit) = e*^*A,(0)e-*^*. (6) 

Here, lO"*^ >j and |0'*' >/ are respectively the wave functions of the initial and final nuclei 
with their corresponding energies Ei and Ej. A denotes the average energy A = {Ei — Ef)/2. 
Ak[0) is the Gamow- Teller transition operator ^4^(0) = J2i'^t{^i)k, k=l,2,3. 
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The time integral form of Mqt in Eqs. (§) and (j^) has been the starting point for the 
methods, which avoid the exphcit calculation of the intermediate nuclear states, e.g. the 



Operator Expansion Method (OEM) [jl5|,|l6[' p^j23|] and the S-matrix approach [^. The 
S-matrix approach requires the derivation of two-body operators from the corresponding 
exchange Feynmann diagrams and the calculation of the nuclear transition by using a given 
initial and final nuclear wave functions. 

For an analytical study of the different approximation schemes, it is useful to rewrite 
the transition operator of the nuclear matrix element MAA{t) in Eq- (iD into a infinite series 
of multiple commutators of the nuclear Hamiltonian H and the Gamow Teller transition 
operator ^^(0) with help of 

n times 

Mt) = £ ^^milH^AMU]- (7) 

n=0 ^• 

If the multiple commutator is calculated without approximation for a nuclear Hamiltonian 
consisting of one- and two- body interactions, we obtain an infinite sum of many-body 
operators. This difficulty may be avoided if some approximation schemes are introduced, e.g. 
the QBA, the renormalized QBA or the approximation schemes of the Operator Expansion 
Method (OEM). We shall discuss this point in the next Section. 

If we integrate over the time variable in Eq. (|^) using the standard procedure of the 
adiabatic switch-off of the interaction as t — oo, insert the complete set of the intermediate 
states |1^ > with eigenenergies between the two axial currents in Eq. (|^) and assume 
that the nuclear states are eigenstates of the nuclear Hamiltonian, we get the well-known 
form of Mgt 

^ n 

^En- Ei + pio + kio En- Ei+ P20 + k2o' 

After the usual approximation pio + fcio ~ Vw + fcio ~ A the form of Mgt in Eq. (H) 
is suitable for the calculation within the commonly used intermediate nucleus approaches 
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(INA) to 2/52z/-decay process like QRPA, RQRPA and shell model methods, which construct 
the spectrum of the intermediate nucleus by diagonalization. 



III. THE APPROXIMATION SCHEMES 

A. Single particle Hamiltonian 

Let restrict our consideration to a single particle Hamiltonian, which in the second quan- 
tization formalism takes the form 

Here, and c^^^ (Cp.m and ) are creation (annihilation) operators of proton and 
neutron, respectively and and e„ are single particle energies of proton and neutron states. 
By using of the Eq. (|^ and the anticommutation relation of particle operators for the time 
dependent axial current Ak{t) we obtain 



^r(i) = E e^^-^^-^")* <P,m,|A(0)|n,m„ > c+„^c„,^^ (10) 

p,mp,n,m,n 

and 

[ir(t/2),ir(-V2)] = o. (11) 

It means that if we consider only the single particle part of the nuclear Hamiltonian, the 
nuclear matrix Mgt is just equal to zero. It is however expected since the 2/32i/-decay is a 
second order process correlated by the residual interaction. Without the residual interaction 
only the two-subsequent beta decay processes are possible, if they are energetically allowed. 
Clearly, without the residual interaction it is not possible to construct the spectrum of the 
intermediate nucleus. 

We note that Mqj^ = comes as a result of the cancelation between the direct and cross 
term of Mgt in Eq. (||). If we use En — Ei+ p2Q + k2o = En — Ef — pio — kio and transform 
Mgt in Eq- ® in the integral representation we obtain 
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Mgt = lim i / < Ot|Afc(0)Afc(-t)e-^*(Pi°+'=io-^) + 

ife(t)ifc(0)e'^*('P^°~'=^°'^)|0+ > dt. (12) 

If we suppose Af^{t) ^ A'lf'it), integrate over time variable in Eq. (|l^) and use the anti- 
commutation relation of the particle operator we find a complete cancelation between both 
terms in the r.h.s. of Eq. (|l^). It shows that the single particle operator of the nuclear 
Hamiltonian plays a less important role in the evaluation of 2/52i/-decay matrix elements. 



This situation has not been noticed in Ref. |27| in which the approximation H //■^•p- was 



also discussed and therefore the authors there came to a different conclusion. 

B. The QRPA and RQRPA Hamiltonians 

The INA approach for the calculation of the Mgt in Eq- (§) consists of two QRPA 
diagonalizations related to the initial and final nuclei. The corresponding initial and final 
QRPA Hamiltonians and take the forms 

W = consU + E ^"^'QImQTm. 

•nii JM 

Hf = const J + ^"^^QT^/QIL. (13) 

mJM 

which are connected with two sets of intermediate nuclear states 

|m, JM, t >= Q^Z'Krpa >^ l^f-JM, f >= Q^j7 Kpa >/ (W) 

generated from initial and final nuclei, respectively. Henceforth we use label "i" for initial 
and "f for the final nuclei. Q"^^'f is the energy of the m-th intermediate state and the 
phonon creation operators Q'jm'''^ is defined as followed: 

pn 

Q^j7 = i:iKn)jB^^PnJM) - Y^^jBipnJM)). (15) 

pn 

^^n)j ^{pn)J forwards- and backwards- variational amplitudes, respectively. 

A~^{pnJM) and B'^{pnJM) are respectively boson creation operators of initial and final 
nuclei. 
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The time dependent axial current Ak{t) can be expanded in the QRPA phonon basis as 



follows [p7 



pnm 



Here, u and v are the BCS occupation amplitudes. 

The expression in Eq. (|TB]) allows us to calculate the transition operator of the time 
dependent nuclear matrix element MAAif) in Eq- If we suppose the equivalence between 
both QRPA Hamiltonians in (O) we obtain 



[A2''^\t/2), A^'^'^-tm] = const, (17) 

which implies that the 2/?2z/-decay nuclear matrix element is equal to zero because the 
2/?2z/-decay operator should be at least a two-body operator changing two neutrons into two 
protons. A generalization of the above discussion to the RQRPA approach is straightforward 
and leads to the same conclusion. It means that the suppression of the 2/?2z/-decay nuclear 
matrix element is connected with the fact that it is a higher order effect in the boson or 
renormalized boson expansion of the nuclear Hamiltonian. We can obtain non-zero results 
only if we go beyond the first order boson or renormalized boson Hamiltonians. 

One can ask why non-zero results are obtained in the INA QRPA and RQRPA calcu- 
lations of Mgt- We believe that the following reasons could be the origin of this problem: 
i) The initial and final QRPA ground states are not orthogonal. Therefore, even for a 
constant transition operator non-zero results could be obtained, ii) The particle number 
non-conservation. We note that even the average particle numbers of protons and neutrons 
for the excited states of the intermediate nucleus differ from the correct ones, iii) There 
is a mismatching between the initial and final QRPA Hamiltonians and as a consequence 
the two sets of intermediate nuclear states generated from initial and final nuclei are not 
orthogonal to each other. In the QRPA or RQRPA calculation of the Mgt one arrives at 
the formula: 
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< 0^rpal^fc(0)|lm^, >< >< 1+ | Afc (0) 1 0+ p„ >j 



rj w > / UI IJLL \ '"V / IIL J. 11 L ,. \ lit- IIL - \ \ / Ul ULl L , . 

M^r= T. '- • (18) 



Here, < 11, 111, > is the overlap factor of the intermediate nuclear states generated from 
initial and final nuclei, given by 

This overlap factor has been considered practically in all QRPA or RQRPA calculations of 
2/?2z/-decay process. In the case in which the two sets of intermediate nuclear states deduced 
from initial and final nuclei are identical Eq. ([I8| ) is just the ortho normal condition for two 
QRPA states. However, this is not the case in a realistic calculation and we shall show 
later that we can hardly expect it within the QBA. In addition we note that the phases of 
the two sets of intermediate states are in principal arbitrary. Therefore, it is necessary to 
identify them e.g. by requiring the diagonal elements of the overlap matrix to be positive 
or by requiring the largest component of the wave function for each state to be positive. 

The equivalence of the two sets of the intermediate nuclear states is connected with the 
equivalence of both QRPA Hamiltonians in Eq. ([T3|). Let discuss this point within the QBA. 

The quasiparticle creation and annihilation operators of the initial (a"*", a) and final h) 
nuclei are connected with the particle creation and annihilation (c"*", c) operators by the BCS- 
transformations. As a consequence there is a unitary transformation between the initial and 
final quasiparticles both for protons and neutrons. In the case of proton quasiparticles it 
takes the form: 



(K 










l( 




ul ] 








( 








Up) 





«p / V Up/ \ -vl ul J \bp 



(20) 



Relation (^) allows us, by using the QBA, to rewrite the boson operators of the initial 
nucleus with the help of the boson operators of the final nucleus: 

A^{pnJM) = UpUnB^{pnJM) - VpVnB{pnJM) (21) 
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We note that \u\ ~ 1 and \v\ ^ 0. Therefore, we shall omit further terms proportional to 
vv. Next, we can rewrite Qjm' with the Q~Ijm^ and QJm follows: 

QIm' = E ioT-^'Q^jZ' + vT'^'Q^l (22) 

where, 

'^J — "'p'^n{^(pn)J^(pn)J ^ {pn)J ^ {pn)J ) i K'^'^) 

pn 

vT""' = E ^P^.(^(:;).C). - KnuKM (24) 

It is now straightforward to rewrite if* with the phonon operators of the final nucleus and 
to perform a comparison with the H-^ . The equivalence of both Hamiltonians requires that 
the following relations are fullfiled: 

J2 nTioT"''oT"''' + p-'^/p-^-'/) ^ ny^s^^^,^ (25) 

mi 

= 0. (26) 

rrii 

We know that in the vicinity of the collapse of the QRPA solution the lowest l"*" state of 
the intermediate nucleus plays an important role in the calculation of Mcx as it is strongly 
influenced by the ground state correlations. Therefore, there is an interest to check the 
validity of the expression (pSj) for this state. We note that the above expressions could be 
used also in the case of the renormalized QBA scheme, if we replace the X and Y amplitudes 
with the renormalized amplitudes X, Y (see Ref. ]10|]) and suppose for the renormalized 
factors the following relation T^\pn)i+ /'^{pn)i+ — 1- The lowest energy of the intermediate 
state Vt!^+ for 2/32z/-decay of ^^Ge obtained directly from the final nucleus by the QRPA and 
RQRPA diagonalization and indirectly with help of expression (^) for three different model 
spaces given in Ref. |2^ is presented in Fig. |1]. By glancing at Fig. 0. we see that close 
to the collapse of the pn-QRPA the equivalence of the initial and final QRPA Hamiltonians 
is getting worse with increasing parameter Qpp. It is more apparent if a larger model space 
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is used. We remark that the accuracy of the pn-RQRPA calculation of the lowest 1^ state 
for a physically acceptable value of gpp is considerably better in respect to the QRPA one. 
It clearly shows that the RQRPA offers a more reliable solution. However, the initial and 
final RQRPA Hamiltonians still remain different and this could be one of the reasons of the 
non-zero results obtained in the INA calculation of Mgt- 

For the sake of completeness we also plot in Fig. § the Mgt of the 2/32i/-decay of ^^Ge 
calculated within the INA pn-QRPA and pn-RQRPA approach for three different model 
spaces. We see that in both cases Mgt demonstrates an increased sensitivity to Qpp with 
enhanced model space. We note also that the pn-RQRPA allows to perform calculations 
far behind the collapse of the pn-QRPA. These calculations have been performed with an 
overlap factor of the initial and final states, which has been approximated as follows: 

<iUiL>^[QlLQ^jM^] = o7'^' (27) 

The advantage of this overlap factor over the one of Eq. (|19|) is that now the results are 
independent of the phases of the quasiparticle states, which are in principal arbitrary [For a 
given quasiparticle eigenenergy E, there are two solutions (u,v) and (-u,-v)]. For the overlap 
factor of Eq. ([T9|) this condition is not fullfiled. This can be easily proved by a numerical 
test. 

It is also of interest to evaluate Mgt by using the overlap functions O^^^ and p™'™'-'' 
and only one of the two RQRPA Hamiltonian, e.g. the W. With the help of Eqs. ([T6|) and 
(p^) one has for the transition amplitudes: 

1 



< i+j|i(o)||o+^, >.= ^E\/^UiK<^-)i + (28) 

V "J pn 



< o+,jii(o)iii+^, >= 7^ E [(^x^;)i +«x>i;;)i)(^r"^' 

V "J pnrrii 



+ («X^W1 + ^>n^Wl)n 'IV^Ul- (29) 

The Mgt calculated with help of only one pn-RQRPA Hamiltonian is drawn in Fig. |^. 
We note a large discrepancy between the results obtained with initial and final pn-RQRPA 
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Hamiltonians and those obtained by the standard INA calculation. In this way one sees 
that the 2/52i/-decay matrix elements are very sensitive to the restrictions of the RQBA 
approximation. From the Fig. ^ it follows also that the INA calculation has been mostly 
influenced by the final pn-RQRPA Hamiltonian. It is an indication that the main problem 
of the INA consist in the non-orthogonality of the initial and final ground states. 

The main message of this Section is that the 2/?2z/-decay is a higher order effect in 
the boson expansion of the nuclear Hamiltonian. Several procedures have been proposed 
which outline the importance of higher order effects and try to take into account some 
higher-order QRPA corrections Maybe, the most perspective way is the fermion-boson 



mapping procedure discussed by M. Sambataro, F. Catara and J. Suhonen p9| within a 
schematic model. In this way it is expected that the correspondence between the initial and 
final nuclear Hamiltonians will be improved and more reliable results could be obtained. 

It is worthwhile to notice that the usual strategy has been first to try to reproduce the 
observed 2/52z/ - decay half times within a given nuclear model in order to gain confidence in 
the calculated 2/?0z/ - decay nuclear matrix elements. However, there is a principal difference 
from the nuclear physics point of view. The 0/321/ - decay is not a higher order effect in the 
boson expansion of the nuclear Hamiltonian and therefore the renormalized QBA scheme 
could be sufficient. It is because the nucleons undergoing beta decays in the nucleus are 
correlated by the exchange of Majorana neutrinos and then the decomposition in Eq. @ 
is irrelevant. We remind the reader that the QRPA and RQRPA Hamiltonians have been 
found successful in describing the single beta decay transition. 

C. The Operator Expansion Method 

The OEM is a nuclear structure method for the 2/52z/-decay, which has the advantage to 
avoid the explicit sum over the intermediate nuclear states. There are two different ways to 
derive the 2/52z/-decay OEM transition operators. In an approach proposed by Ching and 
Ho (OEMl) the expansion of the denominators of Mgt in Eq- (0) Taylor series is 
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used. In a different approach proposed by Simkovic (0EM2) the OEM is derived from 
the integral representation of the nuclear matrix element Mgt in Eqs. (H), @ and (J^). It 
has been found that the 0EM2 offers advantages over OEMl as there are no problems of 
convergence in the power series expansion of the denominator, which has been a subject of 
criticism, J. Engel et al. [0. 

In the OEM the 2/32i^-decay transition operator is described by an infinite series of mul- 
tiple commutators of the nuclear Hamiltonian H and the Gamow-teller transition operator 
A(0) [see Eq. (|^) and (|^]. The OEM is based on two main assumptions [p!^ , p!5| : i) It is 



assumed that the kinetic energy operator can be ignored in the resulting commutators and 
therefore the nuclear Hamiltonian is represented only by the two-body interaction terms, ii) 
Only two-body terms are retained in evaluating each commutator and higher order terms are 
neglected. These two approximations seem reasonable. The omission of the one-body terms 
of the nuclear Hamiltonian is justified by the fact that these terms play a less important role 
in the evaluation of Mgt- It is easy to see it if we consider only the one-body part of the 
nuclear Hamiltonian Hq, then A{t) = e'^^°^A{0)e~^^'^^ as well as the commutator [A(t/2), 
A(-t/2)] are one body operators. However, there are no contributions to 2z/2/3-decay from 
a one body operator, as the 2z/2/5-decay operator should be at least a two-body operator 
changing two neutrons into two protons. In respect to the second approximation from the 
discussion in Section 3 it follows that this approximation goes beyond the QBA or renor- 
malized QBA. In the case of the QBA the commutator [A(t/2), A(-t/2)] is just a constant 
but within the OEM the commutator [A(t/2), A(-t/2)] is a two-body transition operator 
changing two neutrons into two protons. So, it is not true that the OEM is an approximation 
to the QRPA as it was believed before PO . 



Recently, the OEM has been reconsidered |23|. It has been shown that the Coulomb 



interaction plays a decisive role within the OEM. Gmitro and Simkovic |]I6| were first to 
introduce a Coulomb interaction term 

Vc = ^ E 9c{r.,) OJ^, 01^ = \{l + r°)(l + r°), (30) 
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in the OEM formalism. Later it was done by Muto pO[, who however obtained different 
OEM transition operators. We shall discuss this important point and explain the origin of 
these differences. 

The time dependent nuclear matrix element MAA{t) in (|^) can be transformed to the 
following form: 

MAAit) = e^(^-^/)" < 0;|e^^"^[A,(V2), Afc(-t/2)]e-^^"|0+ > . (31) 

Obviously M^^(t) does not depend on a, if the nuclear states and their corresponding 
energies can be considered as the eigenstates and eigenvalues of the nuclear Hamiltonian H. 

H\Ot >= EilOf >, H\0f>=Ef\0j> (32) 

Using the machinery of the OEM for the operator inside the brackets in the r.h.s. in Eq. 
(^) one obtains 

MAA{t) = e^(^^-^-f)" <Oj\VYl e^^^^^"'")^"-" x 

r>+V^^^^(t,r„^,n^(n,m),nr(n,m),^„^)|0+ > . (33) 

Here, 11^ and 11^ are projectors onto spin singlet and triplet states and yo^'^ is the two-body 
OEM potential. V denotes the Principle value integration. It is apparent that M4^(t) is a 
dependent for a ^ 0. It means that if the nuclear Hamiltonian does not contain a Coulomb 
term [|l^,|l^' fl^^ ||19[ or contains a Coulomb interaction in the form of Eq. (§) [|1^,^' ^1 



the derivation of the OEM transition operator is inconsistent. For a = t/2 we obtain the 
OEM transition operators of Muto ||2^ and the formulae of Gmitro and Simkovic could 



be reproduced for a = 0, which differ from each other. 

It is worthwhile to notice that the above results point out some more important aspects. 
The analysis is clearly showing that the mass difference of the initial and final nuclei Ei — Ej 
plays an important role in the calculation of 2z/2/5-decay. 

The above mentioned inconsistency of the OEM could be avoided if one considers an 
effective Coulomb interaction term Vc 
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^c = ^E {Ef-E;)Ol^. (34) 

In this way the one-body terms of the nuclear Hamiltonian are not totally neglected. 

Let's consider the approximated two-body nuclear Hamiltonian H containing central Vcn 
and tensor Vtn interactions in addition to the effective Coulomb interaction (the notation 



of Ref . pO| is used) : 



H^Vc + VcN + Vtn, (35) 



where 



2 . 



gTEir.j) Ulitj) + groir.,) Kitj) ) n^(^j) ], (36) 



Vtn = ^ ( 9TNE{nj) + grNoinj) n^(zj) ) Sij. (37) 

Then, within the OEM approximations the infinite series of the commutators in Eq. (^ 
could be summed using the formulae [|ll,|l^ 



e'^^'Ae-'^^' = \{^ + PAP + cos{2gt){A - PAP) + tsin{2gt)[P, A]), (38) 
^igot^^-igot = A-[0, [O, A]] + cos{gt)[0, [O, A]] + isin{gt)[0, A], (39) 

for P^ = 1 and 0^ = 0. Then by performing the integration in Eq. (^) over t we obtain 
the nuclear matrix element Mgt '■ 

Mgt = < 0^\IvJ2t+t; ( V^^"^'^*(r,,) U^itj) + 

V*"^'^*(r,,) WSl) + V*™(r,,) S,, ) |0+ >, (40) 

where, 



et 



gTE — gsE — ^gxNE + A gTE — gsE + '^gXNE + A 



14 4 4 
■^triplet _ t_ I |_ 



3 A -GgTNO + A GgTNO + A 
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V 



9so ~ 9to + '^9tno + ^ 9so — 9to " '^9tno + ^ 
11 1 2 

tensor _ r_|_ _|_ ^ _|_ 

3 A -'09TNO + A Q9TNO + A 

1 1 



9so — 9to + ^9tno + A (750 — 9to — '^9tno + A 
If we neglect the tensor potential Vtat we have 



(41) 



Mgt =< O^llv Y: r^r^^U'ml^ + f\r)nlJ\Qt > (42) 



fir) ^ 



9TE{r) - 9sE{r) + A 



4 2 

^'^"■^^ A"^?5o(r)-^7To(r) + A ^^^^ 



If we neglect both central and tensor potential V^at we obtain 



Mgt =< 0;|^ E ^,n|;^n ■ c?„|0+ > . (44) 



We note that the OEM transition operators obtained by Muto differ from the above 
ones. This is due to the inconsistent derivation in the framework of the OEMl. The approx- 
imate Hamiltonian considered by Muto does not allow to fix the difference of the masses 
of the initial and final nuclei, which is a basic feature of the 2/32i/-decay matrix element. 
Thus, one obtains different transition operators for different values of the parameter a in 
Eq. This important result is coming from the Heisenberg nature of the axial current 

operators and it can been hardly seen within the OEMl method proposed by Ching and 
Ho |lj]. The importance of the mass difference between the initial and final nuclei becomes 
apparent within the 0EM2 , which contains elements of the field theory. In addition, the 
0EM2 shows explicitly that only the Principal value part of the potential is relevant 

for the nuclear matrix element Mgt- A fact which is not clear within the OEMl. We note 
that the poles of the I^^em potential appear for r„m ~ 1-5 — 2.0 fm. It is because of the 
small energy release for this process (A 2 MeV). It means that only the long range part 
of the nucleon-nucleon interaction plays an important role in the calculation and that the 



yOEM potential is expected to be independent of the chosen type of the nucleon-nucleon 
interaction. It is worthwhile to notice that the consistent OEM-potential in Eqs. (|iO|)-(^^ 
does not vanish even for a zero-range 5-force. 

For the calculation of Mgt in Eq. (|39| ) it is necessary to know the wave functions of the 
initial and final nucleus. The OEM could be combined with the ground state wave functions 
of the QRPA or renormalized RQRPA models 13" IH (OEM+RPA) 



^OEM+RPA^ J2 (^-)h+h'+J+J(^2J +1) 



jp in J 
in' jp' J 

f <Oia II [Cp'Cn']j II J^'ruf >< .rrUflJ^'mi >< J'^rrii \\ [c+c„]j || 0+„ >i 



pnpn 



X < :r|Pr+r2+V«^*^(t, ri2, n^(l, 2), n,^(l, 2), 5i,2) |n, n'; J > (45) 

In the previous OEM+QRPA calculations with the inconsistent OEM potential yOEM ^ ^^iq 
been found not sensitive to the strength of the particle-particle interaction 
g. Recently, the Mg|*^+^-^^ has been calculated for the 2/52i/-decay of ^^Ge within 



the renormalized pn-QRPA with the overlap factor in Eq. ( [TP| ) and a consistent OEM- 
potential in Eqs. (^ID -(^Il) ||23| . Small and large model spaces containing respectively the full 
2 — Ahuo the full — bhuj major oscillator shells have been considered. A strong suppression 
of the results with increasing has been found m. The sensitivity of M^^*^+^^^ to 



Qpp within the physically acceptable region of particle-particle strength has been increased 
considerably with the enhancement of the model space. It is clear that this effect could 
have its origin only in the pn-RQRPA wave functions. We recall that the OEM-potential is 
a two-body operator, which represents a sum over all intermediate nuclear states and it is 
independent of the basis. Then, the instability of the results has to be related to the overlap 
factor and to the fact that the initial and final states are not orthogonal. 

We mentioned already that the overlap factor of Eq. (|T^ does not guarantee the inde- 
pendence of the results in respect to the phases of the BCS states, which are arbitrary. We 
have found that the 2/52i/-decay matrix element is very sensitive to this problem. One can 
see it by calculating M^t'^^^^^^ with different overlap factors, i.e. a delta function overlap 

19 



factor and the overlap factors of Eqs. ([I9D and (|27|). The uncertainty of the results is build 
in into the calculation through the pn-RQRPA wave functions. It is worth to notice that a 
similar problem appears by calculating the closure matrix element, 

Mdos =< Yl -^n^^^n ■ ^m\Ot > • (46) 

We note that Mdos differs only by a constant with the Mqem in the case in which central and 
tensor interactions are neglected. From the above discussion it follows that OEM should be 
not combined with QRPA or RQRPA. Nevertheless, there is an interest for OEM calculation 
with other ground state wave functions especially with shell model ones. 

IV. CONCLUSIONS 

In summary, the integral representation of the nuclear matrix element Mgt in Eq. (^, 
(^) and (|^ has been found useful for a critical analysis of different approximation schemes. 
The single-particle approximation of the nuclear Hamiltonian implies Mqt to be equal to 
zero because of the mutual cancelation of the direct and cross terms. It means that without 
the residual interaction there is no 2/32i/-decay and only the two subsequent beta decays are 
possible, if they are energetically allowed. We note that without the residual interaction no 
excited states could be generated. 

The use of the pn-QRPA and pn-RQRPA nuclear Hamiltonians reveals that the 2(32v- 
decay transition operator is a constant, if one assumes that the initial and final nuclear 
Hamiltonians correspond to each other. However, 2/?2z/-decay operator should be at least 
a two-body operator changing two neutrons into two-protons. Therefore, the inclusion of 
higher order terms of the boson expansion of the nuclear Hamiltonian is necessary for a 
QRPA and RQRPA treatment of the 2/32z/-decay process. It is worth mentioning that there 
is no such requirement in the case of the single beta decay and 2/50z/-decay calculations. 

By using a unitary transformation between the quasiparticles of the initial and final 
uncorrelated BCS ground states the problem of the equivalence of the initial and final pn- 
QRPA (pn-RQRPA) nuclear Hamiltonians has been studied. It has been found that the 
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pn-RQRPA Hamiltonians demonstrate a better mutual agreement like the pn-QRPA ones. 
The mismatching between both pn-QRPA Hamiltonians is large and their correspondence 
is getting worse with increasing model space and particle-particle strength parameter Qpp. 

Previously, it was believed that the OEM is an approximation to the QRPA approach 
2^. However, the 2/32z/-decay transition operator within the OEM is a two-body operator 



and not a constant as within the QRPA or RQRPA. It means that the OEM approximations 



go beyond the QBA and RQBA. We note that the previous derivation pO[ of the 2/32i/-decay 
OEM transition operator was inconsistent because the role of the energy difference between 
the initial and final nuclear ground state in the calculation was overlooked. We present the 
OEM potential derived in a consistent way by considering the effective Coulomb interaction 
term. In this way the one-body terms of the nuclear Hamiltonian are not totally neglected. A 
combination of the OEM transition operator with the pn-RQRPA wave functions reflect the 
instabilities incorporated in the pn-RQRPA treatment of the two-vacua problem. Therefore, 
we suppose that the combination of the OEM-transition operator with other ground state 
way functions, e.g. shell model ones, could be perhaps more predictive for the 2/32i/-decay 
process. 

We remark that the OEM is a special method developed for the 2/?2z/-decay, which could 
perhaps find a wider use for the study of the pion p-wave contribution to the process of 
Double Charge Exchange (DCX) pion on nuclei. The application of the OEM method is 
limited to the problems, where the nucleon-nucleon correlations by a meson exchange plays 
a dominant role. However, this is not the case in single beta decay and neutrinoless double 
beta decay calculations (there is a neutrino correlation of two beta decays in the nucleus). 
For these processes the one-body term of the nuclear Hamiltonian is expected to play a 
crucial role. Therefore, methods constructing explicitly the intermediate nucleus spectrum 
could be more successful in the treatment of these processes. 

The presented studies have shown that the two- vacua problem appearing in the calcula- 
tions of the 2/32i/-decay matrix elements are not safely treated within the QBA or RQBA. 
As a result there are two-sets of intermediate nuclear states generated from initial and final 
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nuclei, which do not correspond to each other. It is because of the particle number non- 
conservation and the violation of the Pauli exclusion principle. Both these problems are 
connected to each other because one can not guarantee the particle number conservation 
without an explicit consideration of the Pauli principle. We note that the RQBA takes into 
account the Pauli principle only in a approximate way. In the QRPA and RQRPA INA cal- 
culation of the 2/30i/-decay matrix element the inaccuracy of the considered approximations 
for solving the two-vacua problem is covered by the introduced overlap factor. In this paper 
we have proposed a less critical way to evaluate it. We believe that by including higher order 
terms of the boson expansion of the nuclear Hamiltonian the equivalence of both nuclear 
Hamiltonians and of the two sets of the intermediate nuclear states will be improved and 
the overlap factor will be closer to a delta function. A schematic study of M. Sambataro et 
al , based on the fermion-boson mapping procedure is rather encouraging and could lead 
to a better agreement of the two-different nuclear Hamiltonians and to more reliable results. 
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FIGURES 

FIG. 1. The energy of the lowest 1"*" state of the intermediate nucleus "^^As calculated directly 
within the pn-QRPA and pn-RQRPA for "^^Se isotope and indirectly with the help of Eq. (25) is 



plotted as function of the particle-particle coupling constant gpp for a 12-level (the full 2 - 4^10; 
major oscillator shells) and a 21-level model space (the full — 5fiuj major oscillator shells). 

FIG. 2. The nuclear matrix element Mgt for the 2z^2/3-decay of "^^Ge calculated within the 
pn-QRPA and pn-RQRPA INA, is plotted as a function of the particle-particle coupling constant 
gpp. The dashed line corresponds to the 9-level model space (the full 3 — Ahuj major oscillator 
shells), the dot-dashed line to the 12-level model space (the full 2 — 4:huj major oscillator shells) 
and the solid line to the 21-level model space (the full — 5huj major oscillator shells). 

FIG. 3. The calculated nuclear matrix element Mgt for the 2z^2/3-decay of "^^Ge pn-RQRPA 
INA plotted as a function of the particle-particle coupling constant gpp for a 12-level and a 21-level 
model space. The solid line corresponds to the standard calculation with the initial and final 
QRPA Hamiltonians. The dashed line is the calculation which considers only the initial QRPA 



Hamiltonian and the overlap functions [see Eqs. (|28| ) and (2£)] and the dot-dashed line the final 
QRPA Hamiltonian and the overlap functions. 
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